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Stability criterion for the surface gravity capillary waves in a flowing two-layered fluid system 
' with viscous dissipation is investigated. It is seen that the dissipative instability of negative energy 

^\ , waves is absent,- contrary to what earlier authors have concluded. Their error is identified to arise 

■ from an erroneous choice of the dissipation law, in which the wave profile velocity is wrongly equated 

' to the particle velocity. Our corrected dissipation law is also shown to restore Galilean mvariance 

^ . to the stability condition of the system. 

(D ■ 

Q: 

^\ , Dissipative instability in a flowing system has been discussed in the literature, and is being sought to be applied to 
explain various phenomena for a few decades. This instability was shown by many authors to be caused by viscosity, 
' ^ ' but seen only in some selected frames of reference in which a particular mode of the system possesses negetaive energy 
p]-pl|. Their result, therefore, implies that the stability of such a system is frame dependent, which obviously violates 
'"^ the condition of Galilean invariance. In the present paper, it is shown that the violation of Galilean invariance arose 
due to a wrong choice of the dissipation law. This question has also been recently addressed by us in an earlier 
paper [T^ , where we calculated the total energy in a magnetohydrodynamic shear flow. It is there shown that the 
existing theories give a rate of entropy production which is not invariant under Galilean transformation. The method 



^ ' to calculate the correct rate of energy dissipation is given by us in the earlier paper in considerable detail. In this 



Oh: 



X 



research note, we analyse the problem from the point of view of Euler's equation of motion and show the validity of 



• our earlier conclusions. 



The system considered is a single interface of discontinuity separating two uniform incompressible fluid media. The 
interface is along the x- direction with the force of gravity acting downwards, i.e., in the negative z- direction. The 
equilibrium pressures, densities and velocities are given by 

■ / \ / \ / \ f Pi, Pi, ui, z < 0, , . 
' with the equlibrium pressure balance condition requiring that pi = P2- 

, We consider small perturbations about the above equilibrium configuration. The vertical displacement ri{x,z,t) of 
a streamline at any point (x, z) can then be written in a form 

l> . 

0> . / . _ ( fi{k,t)exp{ikx)exp{kz}, z<0, , , 

■ '^^^''''^' = \fi{k,t)ejiplikx)eTip{-kz},z>0, ^ > 

' 55 ' that is consistent with an incompressible and irrotational flow. From the linearised equations of hydrodynamics, we 
then obtain the expressions for the velocity and the pressure fluctuations. These expressions are (with i =1,2) 

Qh' Suxi — T {n + kui) fi{k,t) exp (ikx) exp {±{kz)}, (3 a) 



Suzi = I fi(k,t) exp (ikx) exp {±(kz)}, (3 b) 

\ n J 

and / p ■ \ 2 . 

5pi = ± j {n + kui) fj{k,t)exp{ikx)exp{±{kz)}, (3 c) 

with the upper sign designating i = 1 and the lower sign designating i — 2. In Eq.(3), we have assumed the temporal 
dependence to be fj{k,t) ~ exp (mt), with a 'dot' designating a time derivative d/dt. 
At the interface (z = 0), the continuity of pressure fluctuations demand 

6pi{x,t) - Sp2{x,t) = -T—r]{x,t) +g{pi~ p2)ri{x,t), (4) 
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with g being the acceleration due to gravity, and T being the suface tension of the interface. While solving for the 
temporal Fourier amplitudes (i.e., the normal modes) rj{k,n) by substituing (3) in (4), we find that for a non-trivial 
solution to exist, one must have 

(n + A:I7)' =(An)2, (5) 

where, 

(An)^ ^ gk (ai - a^) + k^T' - fc^aiaa (wi - ^2)^ , (6) 

with 



a» =Pi/(Pi+P2),(i = 1,2), (7 a) 

?7 = (aiui + 0121*2) , (7 b) 

""""^ T' - T/ (pi + p2) . (7 c) 

Eqs.(5-7) constitute the familiar dispersion relation for the surface gravity capillary waves that exihibits the Rayleigh- 
Taylor and the Kelvin- Helmholtz instabilities under certain conditions [ |l3[ . 

It is necessary at this stage to introduce the question of dissipative instability. It was argued by several authors 
|^]-|rT| that, on introducing a kinematic viscosity v in any one of the two media (in the lower medium, say), one 
changes the dispersion relation to 

{n + kuf = {Anf + ivaink^ , (8) 
which, for (An)^ > and for a small kinematic viscosity {yaik^U / {AnY l), gives the two roots as 

n± = ~k{U± (An)/fc) + [(An) ± kU\ . (9) 

For U > (An)/k, the (-) root in Eq.(9) is called a negative energy wave that grows in time as 
exp [z/Qfifc^ {kU/ (An) — l) /2] , thus giving dissipative instability of the negative energy wave. Note that, the growth is 
possible for any non-zero but small value of z/ whenever U > {An)/k, while for v exactly equal to zero, the instability 
criterion has no dependence on U and is given by (An)^ < 0, or, 

k'^aia2 {ui — U2)^ > gk {{ai - a2) + k'^T' / g] . (10) 

The above result contains two surprising conclusions. Firstly, the response of the system is not continuous with 
respect to as — > 0, i.e., the stability of the system for an arbitrarily small viscosity is different from that when 
viscosity is exactly zero. Furthermore, the stability with a small but nonzero v appears to be dependent on [/, where 
U is of course dependent on the frame of reference. This means, that the stability of the system depends on the frame 
of reference of the observer. In other words, by moving the observer with a given speed, one can create an instability 
of the negative energy waves,- a result which obviously violates the fundamental law of Galilean invariance. 

The problem can, however, be resolved in the following way. As has been done by most authors, we consider one 
of the fluids, i.e., the lower one to be viscous, while the upper one to be non-viscous. This simplification enables us 
to ignore the complications due to boundary layers, while the essential physics remains unaltered. We note, that by 
substituting in = d/dt and ik = d/dx in the dispersion relation (5) for the non- viscous case, we obtain an equation 
of motion 

-^Tl{x,t) = -{Anfi^{x,t), (11) 



where, D/Dt is the total derivative 
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From Eq.(12), it is seen that, Eq.(ll) describes the force equation of a fluid system moving with a nett velocity U 
with respect to the given frame of reference. Due cognizance should be taken about this fact while calculating the 
viscous force term i/pid'^ /dx^{SvzYlt is to be noted that the term 6vz here denotes the real velocity of the fluid particles 
pertaining to wave propagation [Q, and not the profile velocity dri{x,t)/dt, as the earlier authors have suggested. 
We note that, in this moving fluid system, the particle velocity is calculated as dvz = Dri{x,t) / Dt — i{n + kU)r]{x,t). 

The above argument implies that, the equation of motion in the presence of viscosity must read 



£)2 

— r)(fc, t) = -{Anffjik, t) - uaik'—fiik, t), 



thus giving a dispersion relation 



that yields the two roots 



(n + kU) = {Anf + ivai {n + kU) k^ 



^-Ik^A} 



V a 



1/2 



ivaik^ /2. 



(13) 



(14) 



(15) 



In Eq.(15), the last term on the right hand side predicts a damping for both the wavemodes when (An)^ > v'^a\k'^ /A. 
A growth is, however, possible if and only if 



{Anf - v'^alk^/A < 0, &nd{v^alk^lA - {An)^}^^^ > va^e 12, 



(16) 



thus presenting the same instability criterion (An)^ < 0, as in equation (10) for the non- viscous case. While precluding 
the dissipative instability of negative energy waves, Eqs.(15) and (16) thus suggest that the presence of viscous 
dissipation does not at all alter the stability property of the surface gravity capillary waves. 

The foregoing analysis shows that, the dissipative instability is simply an artifact of an erroneous choice of the 
viscous damping law by the earlier authors. The correct viscosity law, written as in the last term of Eq.(13), 
corresponds to a resistance proportional to v pid"^ {Sv^) / dx'^ = —iupik^ (n + kU) r]{x,t). Replacing this expression 
by —ivpink^rj(x,t) , as was done by the earlier authors, would be equivalent to having a resistance of the form 
a (n, U) V piff^ [Sv^) / dx^ , with a{n, U) = n/ (n + kU). For small values of the kinematic viscosity i/, we can then 
use Eq.(5) to write a(n+,[/) « {kU + An) /An and a(n_,C/) w (— fcJ7 + An) /An for the (+) and the (-) modes, 
respectively. Here, Q;(n+, U) is always positive, but Q;(n_, U) is negative when U > [An)/k. In such a situation, the 
force a(n_, U)v pid^ {Sv^) / dx^ would act as an attractive force, that helps to build up the amplitude of the negative 
energy wave. The erroneous resistance formula, that has been used so far in the literature, thus makes the viscous 
drag force frame dependent and gives a velocity dependent acceleration in selected frames, rather than a deceleration 
in all frames. By substituting unity for a(n, J7), the correct viscous resistance formula however demands that, the 
velocity 5vz be the true velocity of the fluid particles pertaing to wave motion so that, the viscous force becomes frame 
independent. The use of this correct formula leads us to two important results. Firstly, it precludes the possibility 
of dissipative instability in negative energy wave systems. Secondly, it gives a dispersion relation in which the drift 
velocity U appears only as a Doppler shift term in the frequency;- the important consequence of which is that, the 
stability of the system is independent of [/, thus giving the required frame independent universal character to the 
stability condition of the system. 
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